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Abstract. Pressure load Love numbers are presented for 
calculating the mantle deformation induced by the variation of 
the pressure field at the core mantle boundary (CMB). We find 
that the CMB geostrophic pressure fields, derived from "frozen - 
flux ' core surface flow estimates at epochs 1965 and 1975, 
produce a relative radial velocity (RRV) field in the range of 
3mm/dccade with uplift neai the equator and subsidence near 
the poles. The contribution of this mechanism to the change m 
the length of day (l.o.d) is small — about 2.3 x 10 
ms/decade. The contribution to the time variation of the ellipticity 
coefficient J is more important. — about -1.3 x 10 /yr. 

Introduction 

It is generally believed, after Elsasser (1946), Roberts and 
Scott (1965), and Backus (1968), that the secular variation of the 
geomagnetic field (SV) is primarily due to the advection ot pre- 
existing lines of magnetic force of the main field by the fluid 
motion just below the core mantle boundary (CMB). This 
"frozen-flux" mechanism enables us to probe the dynamic 
behavior of the fluid metallic core on the decadal time scale using 
observed SV. However, the inversion of the SV through the 
induction equation for the tangential tlow near the CMB is highly 
underdetermined. One dynamical assumption that helps secure 
uniqueness is that the tlow in the upper reaches of the core is in 
geostrophic balance with the pressure field there (Ball et al., 
T 969; Le Vlouel et al.. 1985. Backus and Le Mouel 1986). 

Pressure variations at the CV1B also cause deformation of the 
solid mantle. Thus, changes in the pattern of flow in the core 
result in deformation of Earth's surface. For world wide geodetic 
networks, this surface deformation is a potential source of signal 
that has not been explored before. If this surface deformation 
were large enough to be observed, and could be separated from 
other effects, it would provide an additional constraint on the 
dynamics of core tlow. The transfer of net angular momentum 
from the core to the mantle via topographic coupling at the CMB 
(Hide, 1969; 1989) is an integrated global signal. In contrast, the 
surface deformation field is only a chart of local signals. Thus, 
pressure variations associated with geostrophic flow, and perhaps 
higher order flows, as well, can not escape influencing the 
deformation field. In this paper we use the geostrophic pressure 
variations derived from core tlow models (Gire et al.. 1990) to 
calculate the radial deformation of Earth's surface. This 
deformation in turn affects Earth rotation by changing the 
moment of inertia. We also examine this two-way interaction by 
calculating the change in the l.o.d and the time variation ot the 
ellipticity coefficient A induced by the deformation field. 
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Pressure load Love number 

The only difference between a mass load and a pressure load is 
that a pressure load has no load mass contribution to the 
perturbation of the gravity potential. The pressure field external 
to the solid earth acts upon the earth as a surface load, thus, the so 
called "pressure load Love numbers" for the external field can be 
easily defined and calculated by slightly modifying Longman's 



(1962) classical mass load theory (see below). The pressure 
perturbation in the interior of the earth, on the other hand, is a 
three dimensional scalar field like a mass perturbation field The 
Green functions for responses to the internal pressure 
perturbation should be similar to the Green functions for the 
responses to the internal mass density perturbation (e.g., Richards 
and Hager, 1984). Practically, it is difficult to estimate the 
pressure perturbation as a function of depth using the SV data, 
but we can estimate the pressure variation near the CMB (see 
below). In this situation, we treat the mantle as a spherical shell 
which is subjected to a pressure load at the CMB, assuming that 
the pressure perturbation within the mantle is only induced by the 
pressure load at the CMB. This problem is similar to the external 
pressure load problem. The scope of this paper is limited to such 
a ' pressure load" problem, so, we use the nomenclature "pressure 
load Love number at the CMB". The nomenclature "internal 
pressure Love number" used by Lefftz and Legros (1992) is more 
proper for describing the three dimensional pressure perturbation 
problem. 

We define pressure load Love numbers at harmonic degree n, 
h n {r), /„(r), and *„(r),as 
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where U n {r), V n (r), and 0 H (r) are the radial, horizontal, and 
gravity perturbation components, respectively (see Longman, 
1962 for details), q is a point force load located at position 
(r 0 ,0,0), g is the outer surface gravity, and p is the average 
density of the earth (note, the point force load q has the 
dimension of a body force). The radial deformation, s r , and the 
gravity potential disturbance, < p r , at radius r are 
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where p is the pressure field, a is the arc length between the 
fixed point ( and the moving point ($', (p ') , and P n is the 
Legendre polynomial. The integral in (2) is over the pressure 
distribution on a unit sphere. 

We see from (1) that the response at r also depends upon the 
radius r 0 where the point load is applied. For the point force 
load q acting on the outer surface of the earth, r = a< these 
pressure load Love numbers for the surface response, r 0 = a , are 
compatible with the mass load Love numbers (Longman, 1962). 
The fluid core comes into effect through the free slip boundary 
conditions at the CMB (e.g., Richards and Hager, 1984). 

Fig. 1 displays different pressure load Love numbers based on 
the seismic earth model PREM (Dziewonski and Anderson, 
1981). In the cases of CMB pressure load, a comparison between 
the CMB response pcc and the surface response pcs shows that 



the mantle acts as a low-pass filter for the surface response pcs, 
with the cut-off frequency at about degree 9. As is well known, 
the frozen flux" approximation only works on a fairly large 
spatial scale (>600 km) (e.g., Roberts and Scott. 1965). The 
implied upper limit in harmonic degree of about 36 is much 
higher than degree 9. Hence, as far as the surface deformation is 
concerned, the "frozen flux" approximation should be accurate. 

Geostrophic pressure field at the CMB 

For a tangentially geostrophic How, the horizontal velocity u and 
the pressure p at the CMB are related through the horizontal 
momentum equation (e.g., Le Moue! et al., 1985) 

2 p b Q cos (H = nxV H /? (3) 

where p is the density of the core, Q. is the angular velocity of 
the earth's rotation, n the radial unit vector, and V H the 
horizontal gradient operator on a unit sphere. Once the How field 
u is known, we can infer p from u by using equation (3) (e.g., 
Gire et al, 1990). Because of the gradient operator V H , a pressure 
field so obtained misses the zero degree term corresponding to a 
uniform dilation or compression along the radius. This is of little 
interest in analyzing the spatial variation of the deformation field. 
For simplicity, we assume that the zero degree terms in the 
pressure field at different epochs are all the same, and thus, set 
the reference for the pressure field such that the degree zero 
terms are zero. 

A number of geostrophic, steady flow models, u, have been 
produced by selectively inverting the SV data collected over the 
period 1840-1990 (e.g., Gire et al, 1986; Voorhies, 1988, 1991; 
Bloxham, 1989: Gire et al., 1990; Jackson et al,, 1993). As noted 
by Hide et al. (1993), most of these steady flow models are 
similar in appearance. Recently, a time-varying flow model has 
been proposed (Voorhies, 1995). A time-varying flow model 
enables calculation of dynamical surface deformation. A simple 
dimensional analysis shows that the contribution of the inertial 
term in the elastodynamic equation is negligible, so, it is 
sufficient to estimate the surface velocity by using the quasi - 
static deformation at successive epochs. In this study, we use the 
primitive, relatively rough steady flow models G6070.3 and 
G8070.3 (Voorhies, 1988). These models were derived by fitting 
the Definitive Geomagnetic Reference Field (DGRF) model from 
1960-1970 and from 1980-1970 respectively. We set the average 
epochs tor the two time intervals to 1965 and 1975 respectively. 
Because of the strong low- pass filter effect, a change in 
smoothness of the model How does not have much impact on the 
surface deformation. Fig. 2 shows the CMB pressure fields 
computed from G6070.3 and G8070.3. 

Results 

Fig. 3 and 4 are the radial deformation fields calculated using 
equation (3). The summation of harmonic degrees starts from 
degree 2; thus, the deformation fields only represent the departure 
from the initial state, presumed to be hydrostatic equilibrium, and 
do not depend upon the choice of the origin of the coordinate. 
Comparing Fig. 2 and Fig. 3, we can see that the small scale 
features in the pressure fields have been filtered out in the surface 
deformation fields, although the pressure fields themselves are 



not very detailed (only up to degree 16). The horizontal 
deformation is, in general, much smaller than the radial 
deformation, and will not be discussed in this paper. 

The relative radial velocity (RRV) fields in Fig. 4 are obtained 
by subtracting the radial deformation at the ” 1 965" epoch from its 
counterpart at the "1975" epoch. A striking feature of the RRV 
fields, both at the CMB and the surface, is that the relative radial 
motion between the "1965" epoch and the "1975" epoch is 
roughly centered around a pole (RRV pole) near ( 60° N, 120°E). 
This "1975” RRV pole is roughly on the opposite side, with 
respect to the North pole, from the 1975 geomagnetic pole at 
(78. 6 U N,70.5° W) (e.g., Rikitake and Honkura, 1985). There is 
compression near the RRV poles and dilation near the equator 
relative to the RRV poles. This structure seems to suggest a 
strong influence of the earth's rotation on the RRV field If, 
instead of the Coriolis force, we put the centrifugal force in (3), 
we would expect a deformation field similar to Fig. 4 but 
centered at the rotation axis. 

The flattening of the figure along the RRV poles increases the 
moment of inertia about the rotation axis and slows down Earth 
rotation. This deformation mechanism differs from that of 
topographic coupling (Hide, 1969, 1989). We calculate the 
change in the l.o.d caused by this two-way interaction mechanism 
using the relation 

<5(1. o.d) = (l.o.d) — — (4) 
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where C m is the mean polar moment of inertia of the mantle, and 
3c^ is the difference in perturbation of the mantle polar moment 
of inertia, c n , between different epochs. Using equation (2), we 
have 
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where wo 0 is the degree 2 zonal harmonic coefficient of the 
CMB pressure field (fully normalized). For the PREM earth 
model. G6070.3 and G8070.3 flow models, and the value C m 
= 7. 2 x I0' 7 kg m~ . we obtain 

_7 

5 (l.o.d) = 2.3 x 10 “ ms / decade (6) 

From (5), we can also calculate the time variation of the 
ellipticity coefficient J ^ using the conservation of the trace of the 
inertia tensor (Rochester and Smylie, 1974): 

/, = -1.3 x l(f" / yr (7) 

Discussion 

The predicted <5(1. o.d) in (6) is two to four orders of 
magnitude smaller than predicted based on the CMB topographic 
coupling mechanism (e. g. Jackson et al, 1993; Hide et al, 1993). 
It is also below the noise level of the observed decadal 
fluctuation in l.o.d (e.g., Lambeck, 1980). On the other hand, the 
predicted J 2 in (7) is quite large. Although the complete 
separation of J 2 from the lumped contribution of ail even zonal 
harmonics is still at issue, there is no question that J 2 dominates 
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the observed satellite nodal drift. If we take the observed J 2 as 
about - 2.5 x ICT 1 1 /yr (Cheng et al, 1989), our prediction would 
provide about 50% of the observed J 2 . This mechanism has been 
overlooked in previous investigations. There are other 
mechanisms, such as postglacial rebound (Peltier, 1985), and 
present-day Antarctic mass changes (James and Ivins, 1995), 
which could predict comparably important ./)■ The uncertainty in 
the present-day secular ice mass changes in Antarctica results in 
predictions of J 2 ranging from -4.1 x 10“' 1 /yr to 1.8x10"" /yr 
(James and Ivins, 1995). Lefftz and Legros { 1992) found that if 
the zonal harmonic component of the CMB topography has an 
amplitude of 15 m and lasts for 5000 years, the viscoelastic 
relaxation of the earth would contribute to J 2 by about 
- 2.4 x 10~' 1 / yr. These previous results, plus our calculations, 
indicate that the observed J 2 is a lumped signal that should not 
be attributed to a single causal mechanism. 

The amplitude of the surface RRV field (Fig. 4) is in the range 
of 3 mm/decade. Although a signal of this magnitude is in 
principle observable with a global VLBI network (e.g., Herring 
1995), it is substantially smaller than other geodynamic signals, 
such as post glacial rebound. In fact, 3 mm/lOyr is even smaller 
than the uncertainties in the estimated deformations induced by 
post glacial rebound (e.g., Mitrovica et al., 1994). Thus, at the 
current sensitivity, it is more likely to be lost in the noise of other 
geophysical processes. But we have to bear in mind that the 
theoretical framework for the dynamo is still incomplete. If 
pressure changes associated with the ageostrophic component of 
the How are large, the associated RRV would also increase. But 
we can be certain from this investigation that no matter how the 
theory develops in the future, no detailed information beyond 
harmonic degree 9 about the CMB pressure variation will be 
contained in potentially observable surface deformation. Finally, 
our analysis of the RRV is based on the steady geostrophic flow 
models at the average epochs " 1 965 '' and T975”. There is no 
obvious evidence to assert that the RRV in Fig. 4 will remain for 
a longer period of time. If Fig. 4 is only a snapshot of a rapidly 
changing movie, it will be even more difficult to measure the 
RRV using geodetic means. 
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Fig. 1. The spectra of radial and gravity load Love numbers h n 
and k n .The spectrum of each type is indicated by a thin line. The 
code p means pressure load, m means mass load, s denotes the 
earth's surface, e denotes the core mantle boundary (CMB). Thus, 
pcs represents the pressure load at the CMB and surface 
response, and so on. The minus sign indicates the direction of the 
load at the surface. 

Fig. 2. The CMB geostrophic pressure fields at the average 
epochs 1965 and 1975. The tangential velocity models used for 
the conversion (equation (3) in the text) are from Voorhies 
(1988). Note, there are no degree zero components in these 
pressure fields. 

Fig. 3. The radial surface deformations at the average epochs 
1965 and 1975. 

Fig. 4. The relative radial velocity (RRV) at the surface and the 
CMB respectively, obtained by subtracting the 1965 radial 
deformation from the 1975 radial deformation (see Fig. 3). The 
white dot marks the location of the geomagnetic pole at the epoch 
1975 (e. g., Rikitake and Honkura, 1985). It roughly shares the 
same great circle with the RRV pole Hanking the geographic 
North poie. 
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